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CHARGE-CHARGE  CORRELATION  FUNCTIONS  IN 
ULTRA  DENSE  PLASMAS 


I.  INTRODUCTION 


Theoretical  investigations  into  the  properties  of  dense  plasmas  have 

benefitted  greatly  from  the  results  of  computer  simulations  of  model 

systems.  The  simplest  of  these  models  is  the  classical  one-component 

plasma  (OCP)  which  is  a  system  of  classical  ions  Imbedded  in  a  uniform 

12  3 

neutralizing  background.  Monte  Carlo  »  and  molecular  dynamics 
simulations  of  the  OCP  have  guided  theoretical  pursuits  and  have  provided 

benchmarks  for  analytic  calculations  of  its  properties,  which  are  now  well 

a  5 

understood.  * 

The  OCP  is,  however,  a  very  simplified  model  of  a  real  plasma  and  is 
applicable  only  to  systems  which  are  so  dense  that  the  electrons  are 
completely  degenerate.  When  the  electrons  are  non-degenerate,  they  can  no 
longer  be  pictured  as  a  uniform,  rigid  background,  and  the  OCP  becomes  an 
inappropriate  model.  Under  these  conditions,  one  oust  treat  both  the 
electrons  and  ions  as  particles.  A  reasonable  model  which  does  this  is  the 
two-component  plasma  (TCP). 

In  the  TCP  the  electrons  and  ions  are  treated  as  classical  (quasi-) 
particles  that  Interact  through  effective  potentials  which  deviate  from 
pure  Coulomblc  behavior  at  short  distances  in  a  way  that  simulates  the 
essential  quantum  diffraction  effects.  A  specific  form  for  such  a 
potential  is  the  one  suggested  by  Deutsch.^  If  a  and  6  are  species  labels, 
and  if 

*«8  5  h/(2’"a8kBT)1/2 

where  y  is  the  reduced  mass  of  this  interacting  pair,  then  Deutsch's 
a8 

potential  is 

Z  Z  e2 

•  (1> 

This  potential  remains  finite  at  the  origin,  and,  therefore,  prevents  the 
collapse  of  the  system.  It  is  expected  to  give  reasonable  results  for  non¬ 
degenerate  plasmas  with  temperatures  above  the  ionization  potential  and 
Manuscript  approved  April  8,  1983. 
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coupling  parameters  less  then  about  two,  or 


where  a  •  3/4im  -  rgJa  0  and  aQ  is  the  Bohr  radius. 

The  TCP  is  more  difficult  to  simulate  than  the  OCP,  but  recently 
7  ft 

Hansen  and  McDonald  *  have  published  results  of  molecular  dynamics  (MD) 

analyses  of  a  fully  ionized  hydrogen  TCP  using  the  effective  potential  in 

Eq.  (1).  Among  the  properties  calculated  from  the  simulation  data  is  the 

dynamic  charge-charge  structure  factor  Sqq(1c,u).  As  the  Fourier  transform 

of  the  charge  density-charge  density  time  correlation  function,  sqq(k»u>) 

reflects  in  its  shape  the  spectrum  of  longitudinal  modes  in  the  plasma. 

The  motivation  of  this  paper  is  to  use  the  simulation  data  of  Ref.  8  as  a 

reference  against  which  various  kinetic  theoretical  calculations  of 

o 

SpQ(k,u)  can  be  tested.  A  similar  study7  has  already  been  made  for  the 
OCP;  the  intention  here  is  to  generalize  that  work  to  two  components. 

The  types  of  theories  of  interest  here  are  microscopic  theories  based 

upon  formally  exact  kinetic  equations  derived  from  projection  operator10”1^ 

13—1  ft 

or  other  techniques.  A  A  general  characteristic  of  such  theories  is  the 

separation  of  the  so-called  "memory"  operator  into  a  static  mean-field  or 

Vlasov  term  and  a  frequency  dependent  collision  term.  Calculations  of 

density  fluctuations  in  dense  plasmas,  using  just  such  a  theory,  were 

performed  several  years  ago  by  Linnebur  and  Duderstadt.^  Unfortunately, 

they  were  unable  to  estimate  the  short-range  part  of  the  direct  correlation 

functions  appearing  in  the  memory  operator  and  had  to  approximate  those 

correlations  by  their  Debye-Huckel  limits.  We  are  now  able  to  calculate 

the  direct  correlation  functions  using  the  potential  in  Eq.  (1)  in  the  so- 

18  19 

called  hypernetted  chain  (HNC)  equation.  ’  This  is  an  approximate 
integral  equation  method  for  calculating  static  correlation  functions  which 
has  proven  to  be  accurate  for  plasmas  in  which  the  coupling  parameter  is  of 
order  unlty^,  and  produces  excellent  agreement  with  the  static  properties 
calculated  by  computer  simulations  in  Ref.  8.  The  proposal  here,  then,  is 
to  generate  the  HNC  direct  correlation  functions  and  use  them  in  the  model 
kinetic  theories  described  in  Ref.  17. 

In  Section  II,  we  define  the  quantities  of  Interest  and  review  the 
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basic  kinetic  theory  concepts  involved.  Calculations  based  upon  three 
different  collision  models  are  presented  in  Section  III,  and  in  Section  IV 
we  present  our  conclusions. 


II.  KINETIC  EQUATION 


For  a  hydrogen  plasma  the  dynamic  charge-charge  structure  factor  is 
defined  by 

2 

SQQ(k'u)  *  f"  “  2sel<k.“>  +  See(k,u>]  ,  (2) 

where  the  partial  density  correlation  functions  S  fl(k,ai)  with  a, 6  ■  i,e 

ap 

are  given  by 


Sa6(k’a,)  “  /d3PdV  ^  /dteiast  ±  <fa(£pt)f gCHcp')> 


(3) 


=  Jd3pdp'  /  dtelu,t  S^Cktjpp')  . 


Here 


f  (itpt)  -  fd3r«'ii,r  f  (rpt)  (4) 

a  1  a 

is  the  Fourier  transform  of  the  phase-space  density 

N 

f  (rpt)  ■  Z  6  (r  -  r  “(t))  6  (p  -  p“(t))  .  (5) 

*  j-1  j  J 

In  Eq.  (3),  N  is  the  number  of  ions  in  the  system  and  brackets  indicate  an 

average  over  the  equilibrium  ensemble.  It  is  not  customary  to  calculate 

S  _(k,u>)  directly,  but  rather  to  find  it  from 
ap 

Sag(k,(D)  »  2  Re  /d3pd3p'  3ag(kz;pp')  »  (6) 

where  the  Laplace  transform  function  $  (kz;pp')  is 


(7) 


Sag(kz;  ppO  -  f“  dteizt  sag(kt*PP')  . 


The  advantage  is  that  the  functions  S  can  be  found  from  the  solution  to 

06  17 

a  coupled  set  of  kinetic  equations  of  the  form 


(z 


~  Sag(kz5PP')  +J^MC 

a  a 


(P) 


I 

Y-e,i 


c  (k)  fd' 
ay  J 


S  a(kz; 
YS 


P  P  ) 


(8) 


~  I  fd3 p"  4  (kz;  PP")  S  (kz;p*p') 


i  S  (k,t“0;*pp')  . 

otp 


In  Eq#  (8),  M  (p)  represents  the  Maxvell-Boltznann  distribution  function 
a 

(normalized  to  unity),  and  the  c  'a  are  the  direct  correlation  functions 

otB 

defined  by  the  Ornstein-Zernike  relations 


where 


Ca6(k)  “  h~«(k)  "  I  C„v<k)  \«<k>  * 


06 


Y-e,i 


ay 


yfi 


(9) 


Vk)  -  Vk)  '  so8 


(10) 


and 


Sa8(k)  "  f  d3Pd3P'  SaB(k,t  “  °5^)  (11> 

is  the  static  partial  structure  factor  for  species  a  and  0.  (Note  that  in 
Eq.  (9)  and  Eq.  (10)  a  factor  of  density  has  been  absorbed  into  the  usual 
definitions).  The  effects  of  collisions  are  contained  in  the  operators 
4  (kz;pp').  Formally  exact  expressions  exist  for  these  operators,  but 

otp 

they  will  not  be  presented  here.  Instead,  we  will  simply  present  results 
based  upon  various  models  for  them. 

Before  preceding  we  note  for  future  convenience  that,  in  analogy  with 
the  single  species  case^,  Eq.  (8)  can  be  solved  formally  to  obtain  the 
Laplace  transformed  functions, 

Sag(kz)  -  /d3p  dV  S(kz;pp')  (12) 
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la  terms  of  the  functions 


J’a6(kz)  "  /dPdP'  dag(k*5PP")  (13) 

where  the  J  * s  are  solutions  to  the  simplified  kinetic  equations 
op 

<2  "  j.a(ks;pp')  -  I  /dp“  *  Y(k*;pp">  jyg(k*;p"p')  * 

a  p  y*e,l  '  TP 

Ma(p)  fi(p-p')  6ag  .  (14) 

Specifically,  the  results  for  a  system  of  electrons  and  Ions  are 

-  Ecfer  [Vkl>  <J«.<k2>  s«‘k>  +  3.i<k2>  si.'k» 

+  Ccel(k)  -  z(Jefi(kz)  cel(k)  +  Jil(kz)  c^k))) 

x  (^(kz)  See(k)  +  Ju(kz)  Sie(k))  ]  (15) 

and 

s".i(k2) '  ithy  tVk2>  <J'u<k2>  su(k)  +  j'u<k2)  see(k>> 

+  (cel(k>  -  ztf^kz)  cie(k)  +  Jle(kz)  c#e(k))) 

x  (5##(kz)  Se#(k)  +  J#1(kz)  Ste(k))]  ,  (16) 
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where  the  function*  f.  (kz)  and  S(kz)  are  defined  by 

a 

E  (kz)  -  1  -  c  (k)  +  z  (5  (kz)  c  (k)  +  }  (kz)  c  (k)),  a  *  B  (17) 
a  ota  aa  aa  ap  fsct 

E(kz)  -  Ee(kz)  Et(kz)  -  (<=el(k>  -  *  (Jge  cel(k)  +  Jel  cu(k))) 

x  (c£e(k)  -  z  (Ju  cie(k)  +  Jle  cee(k)))  .  (18) 

Various  models  for  the  collision  operators  will  generate  different 
approximations  to  the  J  'a.  Using  these  In  Eqs  (15-16),  and  their 

Ctp 

counterparts  within  the  labels  "e"  and  “1"  interchanged  will  give  estimates 
to  the  dynamic  structure  factors.  These  collision  models  will  be  studied 
in  the  next  section. 


III.  THE  COLLISION  MODELS 

In  this  section  we  will  examine  the  charge-charge  structure  factors 
predicted  by  three  different  models  for  $  and  compare  them  to  the  MD 

ap 

results  of  Ref.  8.  All  three  models  will  require  as  input  the  direct 

correlation  functions  c  (k)  which  will  be  obtained  from  the  solution  to 

a8 

the  two-component  HNC  equation  with  the  potential  in  Eq.  (1).  The  first 
and  simplest  model  we  will  investigate  is  the  collionleas  or  generalized 
Vlasov  model  in  which  $  *  0.  The  second  model  employs  a  simple  Fokker- 

OA 

Planck-like  collision  operator  developed  by  Lenard  and  Bernstein  .  The 

21 

final  model  is  one  suggested  by  Duderstadt  and  Akcasu  which  incorporates 

the  exact  high-frequency  behavior  of  *  a  and  models  its  time-dependence  by 

ap 

simple  exponential  decay. 


A.  The  Generalized  Vlasov  Model 


In  this  model  we  simply  neglect  $  in  Gq.  (14)  to  obtain 

aB 

Vt2)  ■  (1  +  *a(k2»  h  • 


where 


,  k.p/m  M  (p) 
y  (k,z)  -  /dp  a  3 


(20) 


z  -  ic.p/m 

a 

is  the  plasma  response  function.  Substituting  Gq.  (19)  into  Gq.  (16) 

yields  the  usual  RPA  results  for  S  (kz)  with  one  important  difference. 

aB 

The  results  presented  here  contain  the  exact  static  correlation  functions 
rather  than  their  Debye-Huckel  limits. 


Taking  these  functions  from  HNC  data,  we  calculated  Sqq(k,  o>  )  for  r  • 
0.5  and  rg  ■  0.4  for  the  same  values  of  q  (“ka)  used  in  the  MD  runs.®  The 
results  are  shown  in  Fig.  1  where  the  standard  RPA  curves  are  presented  for 
comparison  (the  ordinate  has  been  scaled  by  a  factor  of  100).  It  appears 
quite  clear  that  even  at  the  smaller  values  of  the  wavenumber  q,  where  the 
collective  modes  will  begin  to  exhibit  evidence  of  dissipation  not 
available  in  this  collisionless  model,  that  the  generalized  Vlasov  results 
are  in  good  agreement  with  the  MD  data. 

Axchough  this  mean-field  approximation  may  provide  an  adequate  model 

of  Sqq(k,  (u)  for  T  <  0.5  (as  long  as  q  is  >  0.78),  it  does  not  give  a 

good  representation  of  the  charge-charge  spectrum  at  r  m2  unless  q  is 

restricted  to  much  larger  values.  For  appropriate  comparison  with  MD 

results  at  this  value  of  the  coupling,  we  must  introduce  a  reasonable 

approximation  for  $  in  Gq.(14)  and  solve  for  J  . 

aB  aB 


B.  The  Fokker-Planck  Model 


A  simple  model  for 
frequency,  v, 


$  that  Introduces  an  interparticle  collision 

oB  20 
and  allows  an  analytic  solution  to  Gq.(14)  is 


(21) 


*  FP/.  +-*■' »  >  o,  3  /  3.S  +  \»/+  *'  \ 

•a#  <kz;pp  >  “  1  T  5agVag  —  *  ("7  +  ^  )  «(P  "  P  )- 

3p  3p  a 


Substituting  this  approximate  collision  term  into  Eq.  (14)  and  Fourier 
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transforming  in  momentum  allows  one  to  find 


VkZ)  “  6ag  *a(kz)  * 


(22) 


where 


K  (kz)  -  -  (igm  u  /k2)  I  (Q2,  s  -1) 
a  a  eta  a  a 


(23) 


aDd  2  _2(Ql  +  8a)  Qa  Qa  -x  (Qa  +  Sa) 

X  «T,»  )  -  Q  a  °  e  a  /  °  dx  e  x  x  a  a 

a  a  ci  o 


(24) 


2  2  2 

In  the  above,  Q  ■  k  / gm  v  and  s  *  -  iz/v 

a  a  aa  a  a a 


Having  Introduced  two  collision  frequencies,  vee  and  v^,  into  our 
solution  for  SqQ(k,  u>  )  to  allow  for  collisional  damping,  we  must  now 
select  them  with  some  care.  It  is  well  known  that  the  ordinary  Spitzer- 
like  binary  collision  times  are  inappropriate  in  moderately  to  strongly 
coupled  plasmas.  We  can,  however,  utilize  transport  coefficients  found 
directly  from  the  molecular  dynamics  simulations  to  obtain  estimates  of  the 
required  collision  frequencies.  In  particular  using  the  self-dif fusion 

g 

coefficients  given  by  Hansen  and  McDonald  ,  these  frequencies  are  seen  to 
be 

MD 

vee  _  1  _  1 

^pe  m  gD  3rD 
e  e  e 


.  and 


MD 

Vii 


£>1/2 

®i 


(25) 
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Collisions  are  found  to  be  of  negligible  importance  for  SpQ(k,(a)  at 

r  -  0.5.  They  are,  however,  necessary  when  describing  S^^k.ai)  at  larger 

values  of  coupling*  In  Fig.  2  we  compare  the  collisionless  generalized 

Vlasov  model  calculations  (using  "correct”  HNC-evaluated  direct  correlation 

MD 

functions)  with  the  results  employing  v  from  Eq.  (25)  in  the  Fokker- 

cux 

Planck  solution  for  the  case  r  *  2,  rg  ■  1.  From  the  figure  it  is  clear 
that  the  model  including  collisions  provides  a  much  superior  fit  to  the  MD 
data  than  the  collisionless  model.  For  values  of  q  smaller  than  0.78,  the 
discrepancy  is  still  larger,  indicating  the  failure  of  mean-field  (RPA) 
descriptions  of  strongly  coupled  systems  in  the  collective  regime. 

Even  though  the  FP  model  considerably  improves  the  shape  of  the 

spectrum,  electron  peak  Intensity  for  the  FP  curve,  which  is  a  sensitive 

function  of  v  ,  falls  below  the  MD  data  points  at  peak.  In  addition  the 
act 

peak  position  is  lower  in  frequency  than  is  shown  by  the  simulation.  To 
correct  this  deficiency  we  must  consider  a  different  collision  model. 

C.  Duderstadt  Akcasu-Llnnebur  (DAL)  Model 

Neither  of  the  first  two  models  discussed  above  satisfy  the  third  non¬ 
vanishing  frequency  moment  of  Sp^k.w).  This  sum  rule  will  however,  be 
automatically  satisfied  by  any  model  collision  operator  which  Incorporates 

the  known  high-frequency  behavior  of  4  .  One  such  model  is  that 

21 

suggested  by  Duderstadt  and  Akcasu  ,  and  later  applied  to  weakly  coupled 

two-component  plasmas  by  Linnebur  and  Duderstadt^,  in  which  the  time 

dependence  of  4  .  is  modeled  by  two  relaxation  times.  We  define  the  DAL 
ap 

collision  operator  as 


D  (0) 

DAL,.  ++,.  .  a  3  ,3  .8  ^ 

•«b  (kz:pp  > ' 1  s.s  ~rr:  —  +  r p>  5  <p'p  > 


z+ia^k)  3p  3p  a 


(26) 


2 - M  (P)  P  •  A  (k)  •  p'  , 

c  o  ap 


z+ia  a(k) 
ap 


where 


V0)  ■  1?  /  df  Vr>  • 


(28) 


and 


08 


(k)  ‘  srs-  c<,s(k>  -  £r  f df  *„8<r>  elk-t  ??',„6<t)  • 

a  8  a  8 


(29) 


s  d 

The  decay  constants  a  _(k)  and  a  .(k)  are  chosen  to  obtain  the  correct 

ad  aB 

short  and  long  time,  as  well  as  the  known  large  and  small  wave-number 

limits  of  S  .(k,<d)« 

<*B 

Using  Eq.  (26)  in  Eq.  (14)  and  once  again  Fourier  transforming  in 
momentum  leads  to 
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where  I  is  defined  as  in  Eq.  (23),  but  with  v  replaced  by  w  (kz) 
a  aa  aa 

which  is  given  by 

8  Dn(0) 

w  (kz)  -  -2 - 2 -  . 


aa 


“a  z+iaS  (k) 
aa 


(32) 


The  remaining  functions  appearing  In  Eq.  (31)  are 


yo8(kZ)  S 


k  .  A^k)  .  k 


(33) 
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and 


A  =  zl  -1 
a  a 


(34) 


The  form  of  the  damping  functions  are  chosen  to  be 


*“ed<k>  11  +  'k/^d)2l  • 


(35) 


where  the  unknown  parameters  a  (0)  and  k 

°®21  22 

and  u  constraints  on  S  (k,w).  * 

aB 


s,d 

aS 


are  determined  from  known  k 


With  these  restrictions  we  are 


required  to  choose  only  two  other  collision  frequencies  to  complete  the 
description  of  the  DAL  model.  Those  coefficients  are  the  electron-ion  and 
ion-electron  collision  frequencies,  v  ^  and  ,  characterising  momentum 
transfer  between  particles. 


We  can  find  v  ^  from  the  coefficient  of  electrical  conductivity 


provided  by  the  MD  data. 


MD 


U) 

p 

*e 
4v  a 


1 

4  wo* 


(36) 


As  an  alternative  to  using  the  simulation  derived  conductivity  to  find 
we  can  employ  a  quantum  kinetic  theory  treatment  of  plasma  time-correlation 
functions which  leads  directly  to  the  calculation  of  the  electrical 
conductivity. This  treatment  can  be  thought  of  as  a  quantum 
generalization  of  the  fully  renormalized  kinetic  theory  of  Mazenko1^  in 
which  the  disconnected  approximation  (DA)  is  used  to  renormalize  the 
potential  terms  in  the  collision  operator.  The  application  of  this 
approximation  to  strongly  coupled  plasmas  has  been  discussed  in  Ref.  22  and 
elsewhere. 2^,26  ^ore  recently  the  theory  has  led  to  numerical  evaluations 


DA 


of  o  (or  v  . ,  which  we  label  v  . 

Cl  Cl 

here.  '  We  can  estimate  v  by  taking 


)  for  the  cases  under  consideration 


DA 


MD 

'*ee 


DA 


ee 


MD  vei 
rfai 


The  collision  frequencies  are  collected  in  Table  I  for  the  cases  used 
in  the  simulation  runs.  The  results  employing  the  DAL  model  with  MD  and  DA 
collision  frequencies  for  a  r  ■  2,  r  ■  1  TCP  are  given  in  Fig.  3.  At  low 
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k,  it  Is  appareaC  chat  the  fits  are  superior  to  the  F?  model  In  peak 
position  as  well  as  the  Intensity  of  the  resonance.  At  small  wavelengths 
the  generalized  Vlasov  solution  Is  nearly  Identical  with  those  of  the 
collision  models.  Indicating  the  lack  of  Importance  of  collisions  In  this 
regime. 

g 

In  their  paper,  Hansen  and  McDonald  developed  a  memory  function 
theory  (mft)  based  on  the  work  of  Abramo,  et  al.  to  describe  the  charge- 
fluctuation  spectra  obtained  in  the  MD  simulations.  This  theory  builds  a 
hierarchy  of  memory  functions  for  linear  combinations  of  momentum- 
integrated  microscopic  particle  densities  defined  in  terms  of  known  sum 
rules.  The  hierarchy  is  truncated  at  the  highest  order  in  which  the  sum 
rules  can  be  exactly  calculated  using  two-particle  radial  distribution 
functions  (the  fourth  order  sum  rule).  The  highest  order  memory  functions 
are  damped  in  time  by  an  exponential  or  Gaussian  approximation  employing 
relaxation  times  defined,  again,  by  the  known  sum  rules. 

The  advantage  of  “mft"  lies  in  the  fact  that  only  two  relaxation  times 
are  required  and  that  these  constants  are  determined  within  the  framework 
sum  rule  calculations.  The  DAL  model  on  the  other  hand,  requires  the 
determination  of  four  independent  relaxation  times  by  introducing  four 
transport  coefficients.  However,  since  these  time  constants  interpolate 
between  hydrodynamic  and  Vlasov  behavior,  a  knowledge  of  the  needed 
transport  coefficients  (which  can  be  accurately  found  from  MD  simulations 
or  quantum  kinetic  theory)  ensures  correct  small  k,  small  u>  limits 
of  SqQ(k,(n).  The  charge-fluctuation  spectra  predicted  by  the  DAL  model  and 
mft  are  very  similar  for  all  of  the  cases  under  consideration  here  (see 

Figs.  6  and  8  of  Ref.  8).  Results  of  S  (k,  u>  *  0)  for  these  models  are 

act 

collected  in  Table  II. 

As  noted  above,  one  may  expect  differences  to  occur  at  longer 

wavelengths  where  the  DAL  model  has  used  transport  coefficients  to  ensure 

correct  hydrodynamic  behavior.  In  Fig.  4  we  present  a  comparison  of  mft 

with  the  DAL  model  using  values  of  the  collision  frequencies  obtained  form 

MD  transport  coefficients  and  the  quantum  theoretical  treatment  with  the 

disconnected  approximation  at  q  ■  0.307.  The  mft  results  were  taken  from 

MD 

Figs.  9  and  10  of  Ref.  8.  The  DAL  model  using  v  predicts  a  plasma  peak 
of  greater  intensity  and  lower  frequency  than  the  mft  calculations.  The 


16 
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•  t 


DA 

DAL  model  using  the  larger  values  of  v  predicts  a  smaller  intensity  and  a 
peak  position  closer  to  mft.  Differences  between  these  theoretical 
approaches  can  be  expected  to  increase  as  q  is  made  smaller  still. 
Unfortunately  there  are  no  MD  results  at  these  wavenumbers  to  indicate 
which  is  the  preferable  theoretical  prediction. 

IV.  DISCUSSION 


The  collision  approximations  we  have  investigated  for  describing  the 
TCP  all  utilize  equilibrium  correlation  data  as  input  to  the  solutions. 

O 

This  data  has  been  shown  to  be  accurate  and  easily  obtained  from  the 
solution  of  the  hypernetted  chain  integral  equations  using  the  Deutsch 
potential  as  the  effective  Interparticle  interaction  potential.  We  have 
used  this  data  in  all  models  considered  with  the  exception  of  the  usual 
Vlasov  approximation. 

Q 

As  has  already  been  shown  ,  the  standard  Vlasov  equation,  which 
contains  neither  collisional  dynamics  nor  correct  initial  time  correlation 
information,  does  not  reproduce  any  of  the  MD  spectra  accurately.  The 
generalized  Vlasov  equation,  in  which  the  Fourier  transform  of  the 
potential  in  the  mean-field  term  is  replaced  by  the  exact  direct 
correlation  function,  provides  a  reasonable  fit  to  the  data  at  r  ■  0.5  for 
all  rs  and  q  values  simulated.  This  approximation  is  also  adequate  for 
large  wavenumbers  at  higher  values  of  the  coupling  parameter. 

At  the  higher  value  of  r  considered  here,  it  is  necessary  to  Include 
collisions  in  order  to  damp  the  strongly  spiking  plasmon  peak.  The  Fokker- 
Planck-llke  Lenar d-Bernstein  model  approximates  these  effects  by 
introducing  two  self-collision  times  into  the  solution  of  Sqq(k,u)).  Even 
though  these  constants  are  well-known,  the  resulting  spectra  are  only 
qualitatively  simular  to  the  MD  data  at  smaller  values  of  q.  As  we  have 
noted,  the  FP  term  does  not  satisfy  the  fourth  frequency  moment  sum  rule. 

In  addition  the  collisional  invariant  of  momentum  is  not  satisfied  by  this 
model. ^ 

The  DAL  model  does  incorporate  exact  sum  rules  up  to  fourth  order  and 
does  satisfy  conservation  of  momentum,  leading  to  much  improved  results  at 
low  q  in  the  case  of  large  coupling. The  implementation  of  several 
relaxation  times  allows  the  collision  dynamics  to  be  Interpolated  from 
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known  hydrodynamic  (k  +  0,  u  +  0)  forms  of  S  .(k,u>)  Co  their  Vlasov  and 

ctp 

free  particle  forms*  This  approximate  form  of  the  collision  term  depends 
on  the  accuracy  of  the  acquired  hydrodynamic  coefficients,  but  is  seen  here 
to  give  a  good  description  of  the  MD  charge-fluctuation  spectra.  This 
description  la  very  similar  to  the  mft  calculations  except  at  low 
wavenumber  where  the  interpolative  DAL  model  should  represent  the  spectra 
resonances  as  accurately  as  the  transport  coefficients  provided.  For  all 
cases  r  <  2  and  q  >  0.78,  the  distinctions  are  negligible. 

The  primary  limitation  in  this  process  is  the  accurate  calculation  of 
the  radial  distribution  functions  from  the  HNC  procedure.  This  implies 
that  the  best  solutions  will  arise  from  the  most  inclusive  effective 
interionic  potential.  The  Deutsch  potential,  Eq.  (1),  is  not  valid 
for  r  greater  than  about  two,  but  the  HNC  scheme  should  give  reasonable 
distribution  functions  for  a  more  general  effective  potential.  From  these 
functions  relaxation  times  can  be  derived  and  used  in  the  DAL  collision 
term  to  provide  accurate  predictions  of  dynamic  correlation  functions  in 
more  strongly  coupled  two-componenc  plasmas  over  a  wide  range  of  k  and  u>. 
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